The shape of the A baryon in a covariant spectator quark model 
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Using a covariant spectator quark model that describes the recent lattice QCD data for the A 
electromagnetic form factors and all available experimental data on 7A'^ — > A transitions, we analyze 
the charge and magnetic dipole distributions of the A baryon and discuss its shape. We conclude 
that the quadrupole moment of the A is a good indicator of the deformation and that the A'*" charge 
distribution has an oblate shape. We also calculate transverse moments and find that they do not 
lead to unambiguous conclusions about the underlying shape. 



I. INTRODUCTION 

The determination of the shape of a baryon is an inter- 
esting problem that addresses complex technical issues, 
both from the experimental and theoretical point of view 
P, Recently, for example, there has been particu- 
lar interest (and controversy) about possible deviations 
from a spherically symmetric shape of the nucleon. In 
this work we will address another interesting question, 
namely whether the A baryon has a spherical shape or 
not, and maybe more importantly, what we mean by 
that. 

Information on a baryon's shape is encoded in its elec- 
tromagnetic form factors, which can be measured, at 
least in principle. The electromagnetic form factors de- 
scribe how a particle interacts with a photon, and the 
number of independent form factors depends on the par- 
ticle's spin. A particle with spin-1/2, such as the nucleon, 
is characterized by only two form factors, namely the 
electric charge Geo and magnetic dipole Gmi, whereas 
a particle with spin-3/2, such as the A, has altogether 
four form factors, namely, in addition to Geo and Gmi, 
also electric quadrupole Ge2 and magnetic octupole Gms 
form factors. 

Higher-order form factors and moments can provide a 
measure of the deformation of an extended particle. For 
instance, non- vanishing values of GE2iO) and Gm3(0) in- 
dicate a deviation of the charge and magnetic dipole dis- 
tributions from the spherically symmetric form [3[ . How- 
ever, this simple rule does not apply to pointlike parti- 
cles: Ref. Q shows that structureless spin-3/2 particles 
with charge eA (specified in units of the proton charge) 
possess the "natural" quadrupole and octupole moments 



besides 



G£;2(0) = -3eA, GM3(0) = -eA, (1.1) 



Gi;o(0)=eA, GAn(0)=3eA. (1.2) 



Therefore, the inference of deformation based on non-zero 
quadrupole and octupole moments has to be restricted to 
spatially extended particles. 

It is important to distinguish between "spectroscopic" 
and "intrinsic" moments, the former being the observable 
values of the corresponding electromagnetic form factor 
at zero momentum transfer, whereas the latter refer to 
quantities calculated from charge or magnetic density dis- 
tributions, which are not directly observable [HQl- It is 
only in the nonrelativistic limit that the spatial distri- 
butions of the charge and magnetic densities are related 
to the electromagnetic form factors through a Fourier 
transform, such that information about the shape of the 
density distributions can be accessed through form fac- 
tor data: in the Breit frame, the charge distribution is 
the Fourier transform of Geo, and the magnetic dipole 
density is the Fourier transform of Gmi [E, Q- 

However, this interpretation has its limitations. In a 
relativistic description, the spatial distribution of charge 
or magnetic densities depends on the reference frame. 
When an absorbed photon imparts only a small momen- 
tum to the baryon, this frame dependence can be ignored, 
but in the general relativistic case there is no direct rela- 
tion between form factors and coordinate-space densities. 

To address this difficulty, other concepts and measures 
of shape and deformation were proposed, such as trans- 
verse densities and moments [7HlOl|. and spin-dependent 
density distributions [Tlj . Transverse densities are calcu- 
lated in the infinite momentum frame and give informa- 
tion about density distributions in the transverse plane 
only. This gives valid information in the relativistic case, 
but it is not clear how much one can deduce about the 
intrinsic deformation of a density distribution in a parti- 
cle's rest frame. 

A spin-1/2 particle like the nucleon does not possess 
an electric quadrupole or magnetic octupole moment, 
and therefore it cannot be deformed if finite higher spec- 
troscopic moments are used as criteria for deformation 
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[H, other hand, its density distribution may 

not be spherically symmetric j_giving rise to non-vanishing 
intrinsic higher moments P, 

The situation is somewhat simpler in the case of the 
A, the lightest baryon with spin-3/2 and first candidate 
for a non- vanishing electric quadrupole moment, because 
its intrinsic and spectroscopic quadrupole moments dif- 
fer only by a constant factor [l|, [l^l- Thus, there is a 
strong motivation to estimate the A electric quadrupole 
moment. Previous studies on the deformation of the A 
can be found in Refs. 

In this work we used relativistic quark-diquark model 
wave functions of the A baryon that were constructed 
within the covariant spectator theory [l^ to calculate 
charge and magnetic density distributions, both in mo- 
mentum and coordinate space. Thus we have direct 
information about their shape. In particular, we com- 
pared two models, one which includes only iS-waves and 
is therefore spherically symmetric, and another that in- 
cludes D-wave components that induce a small deforma- 
tion. We determined then how this shape information 
manifests itself in the higher moments, as well as in the 
corresponding transverse moments. 

With a spherically symmetric spatial wave function, 
where the quark and diquark are in a relative S'-wave 
[201, we obtained 



G£2(0) = 0, Gm3(0)-0. 



(1.3) 



This result does not depend on any specific model pa- 
rameters, but holds in general as long as only 5-waves 
are present. 

Once higher orbital angular momentum components 
are included in the quark-diquark wave function, non- 
vanishing values of Ge2{0) and G'm3(0) are generated 
[21] . [22! . This indicates that it is possible to relate devia- 
tions from spherical symmetry in the charge (dipole mo- 
ment) distribution to the value of the electric quadrupole 
(magnetic octupole) moment, similar to what was found 
in the nonrelativistic limit. The calculated values of 
the transverse quadrupole and octupole moments on the 
other hand do not seem to give any clear indication on 
the deformation of the density distributions. 

This paper is organized as follows: The covariant spec- 
tator quark model is introduced in Sec. [ill In Sec. IIII[ 
we relate the electromagnetic form factors with the ex- 
perimentally accessible polarized helicity amplitudes. In 
Sec. IIVI we discuss the deformation of the A as deter- 
mined from transverse densities, and in Sec.|V]we present 
the results obtained from the usual three-dimensional 
densities. In Sec. IVIIwe draw our conclusions. 



II. SPECTATOR QUARK MODEL 

We apply a quark model obtained in the covariant 
spectator formalism [l^, E^, l23l - [3^ . and parameterized 
to describe the A baryon, as discussed in detail in 



Refs. |21|, 122, 133H36[. The A wave function is a mix- 
ture of an S state (L = 0) and two D states (L = 2, 
coupled to core spin 1/2 and 3/2) for the quark-diquark 
system [13, HH , of the general form 



(2.1) 



In this equation, a is the admixture coefhcient of the 
D3 state (quark core with spin-3/2) and b the admixture 
coefficient of the Dl state (quark core with spin-1/2). 
The momentum and spin indices arc suppressed for sim- 
plicity. The S- and D-state wave function components 
are written in terms of spin, orbital angular momentum, 
and isospin operators, multiplied by scalar functions, tps, 
ipD3, and V'-Di- In our covariant spectator model, the 
diquark four-momentum is on-mass-shell, and therefore 
these scalar functions depend only on the square of the 
quark four-momentum, (P — fc)^. where P and k denote 
the A and the diquark four-momentum, respectively. 

Assuming each of the A wave function components in 
(j2.ip to be normalized to 1, the overall normalization con- 
stant becomes N = l/Vl + + 6^. This specific form of 
the wave function was introduced in Refs. [3J,[3^ and two 
different parameterizations were studied in Refs. [2ll.[23|. 
Here we will use the model of Ref . [13 , because it gives 
a more consistent description of the valence quark con- 
tribution of the A to the — >■ A reaction, both in 
the physical region and in the regimes accessible in lat- 
tice QCD calculations with heavy pions. In this model, 
the two ZJ-state probabilities arc both about 0.89% (a = 
0.08556 and b — 0.08572). For more details on the model 
we refer to Refs. 0,113. 

The internal structure of the constituent quarks is 
described in terms of quark electromagnetic form fac- 
tors, parametrized through a vector meson dominance 
mechanism and included in an effective quark current jj^ 



[121 . 130l . I35l . |36| . Employing the wave function (|2.ip and 
the quark current j^, in the covariant spectator formal- 
ism we write the electromagnetic current [l^ [20I [sF 



as 



\ J k 



(2.2) 



Ua{P+,s') 



2Ma 



where P+ {P-} represents the final (initial) four- 
momentum, q = P+ — P- is the transferred momentum, 
Q2 = _qr2^ Ma is the mass of the A, and As the diquark 
polarizations. For the covariant integration over the on- 
mass-shell diquark momentum k we use the abbreviation 



d^k 



(2.3) 



with 



s + k , where mg is a model parameter 
that corresponds to a mean value of the spectator diquark 
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mass [12|, |27[ . The asymptotic states Mq are the Rarita- 
Schwinger vector states [37|. Throughout this paper we 
foUow the convention used in our previous work that the 
diquark polarization indices, Ag, on the wave functions 
are suppressed. 

The multipole A form factors can be written as hnear 
combinations of i^*, i = 1, . . . , 4 [H [M |3|. For 
= 0, to first order in the admixture coefficients a and 
b, one finds [HI 

GEoiO) = N^eA 
Ga/i(0) = N^ca + ka) 
Ge2(0) = 3iaN^)eAl'D3 

Gm3(0) = (eA + KA)iV2[aI^3 + 2&I^i], (2.4) 

where 

eA-i(l+T3), nA = l{K+ + K_f,)^, 

with = diag(3, 1, — 1, — 3), and Mjv is the nucleon 
mass. The factors I'j^^ and I^g are defined in terms 
of the overlap integrals between the initial iS-state and 
the final I?-state, as 

1^3 = lim - f b{k,q,P+)i'D3{P+,k)i;s{P-,k) 
T Jk 

T'di = lim - / h{k,q,P+)i'Di{P+,k)i's{P-,k), 

T Jk 

with T = xatt- The function b{k,q,P+), whose detailed 

form is given in Ref. [s^], reduces to k^y2o(k) in the limit 
— >■ 0, where Y2o{z) is the familiar spherical harmonic. 
The model described above was applied in Ref. [22j 
to calculate the A electromagnetic form factors, and its 
results were compared successfully to the recent lattice 
QCD simulations of Refs. [13, SO]- Although the model 
is still incomplete because important degrees of freedom, 
such as meson (pion in particular) cloud effects, are not 
included, it agrees well with the lattice QCD data for 
Geo and Gmi [13, l4^ a nd is also consistent with the 
unquenched Ge2 data [l3|- This success can be due to an 
effective suppression of pion cloud effects in the elastic^ 
7A ^ A reaction, in contrast to the 7iV A transition 
where the opening of the TriV channel is crucial j33[^35j . 
It is also possible that effective pion cloud effects are 
already included adequately through the vector meson 
dominance mechanism which models the effective quark 
current j^. 



A similar effect can be seen in the nucleon elastic form factors. 
In some models, the pion cloud contributions are around 10% 
[25I l26l | . Even in models where pion cloud contributions to the 
nucleon magnetic moments are significant (fs 40%) [53 . the dif- 
ference between the results with the pion cloud and the results 
when the pion cloud effect is removed differ by only about 5% 

mill. 



III. FORM FACTORS AND HELICITY 
AMPLITUDES 

The electromagnetic form factors of a baryon are in- 
variant functions of Q^. They are independent of the 
reference frame and of the initial or final polarization of 
the baryon. Note that these functions are not directly 
measured in an experiment. What can be measured are 
cross sections in a particular frame, from which helic- 
ity transition amplitudes between two different or equal 
polarization states of the baryon can be deduced. 

In an elastic process like 7 A — !• A, there are only 3 in- 
dependent components of the current, as a consequence 
of current conservation. These components can be cho- 
sen to be J°, and J^, or, alternatively, J°, J+ and 
J-, with J"^ = T^iJ'' ±iJ''^)- Note that is asso- 
ciated with the photon polarizations A = ± that involve 
a change of the baryon polarization (±1), and with 
A = where the baryon polarization is conserved. 

The transition amplitude for spin projections s and s' 

is 

J\s',s) = -u^{P+,s') [0"^'^(e^)^] up{P^,s), (3.1) 

where the operator O"^^ is implicitly defined through 
Eq. (|2.2p . and (e^)"^ are the photon polarization vectors, 
with (e^)0 = (1, 0, 0, 0) and (e^)± = ±^(0, 1, ±i, 0). 

We will work in the Breit frame, where the photon 
four-momentum is q = (0,0,0,(5), with Q = ^/Q^, the 
photon three-momentum q points along the positive z- 
direction, and the initial and final total momenta are 
P± = (AfA\/rTT,0,0,±iQ). 

The spin non-flip components of the current (s' = s) 
are 

J° (s, s) = Geo{Q^) - Ps{s)tGe2{Q^), (3.2) 

for .s = ±i,±f , with /,(±f) = 1 and /,(±i) = -1. 

One can combine the two independent amplitudes as 
a symmetric combination of matrix elements 

4-l[j"{+h+i)+j'{+m)]=GEo{Q% 

(3.3) 

and an asymmetric combination 

J'a - \ [J' i+l+l) J" i+h H)] = -IrGE^iQ'). 

(3.4) 

In the limit — > 0, the first equation yields Geo{0), 
whereas Ge2{0) cannot be obtained directly from the 
amplitudes at = because t goes to zero. 

Similarly, the magnetic form factors are obtained from 
the spin-flip current matrix elements J^(s',s). Again, 



^ In Rcfs. [Toll the normalization is Ua{P, s)u°' (P, s) = —2M/^ for 
P = (Ma,0). Here wc use Mq (P, s)«" (P, s) = -1. 
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there are only two independent amplitudes related with 
Gmi and Cms for ^ (see Refs. [l^ for details). 

In a nonrelativistic formalism, the baryon's shape, and 
in particular any possible deformation — deviation from a 
spherically symmetric form — would depend on the spin 
projection along the z-axis. One can then define an elec- 
tric charge distribution, p£;(r, s), associated with each 
spin projection (±| or ±i) and define an intrinsic elec- 
tric quadrupole momentum [l[ as 



Qa(s) = J d''rpE{r,s)r^Y2oir). 



(3.5) 



When Qa{s) 0, its sign indicates whether the system 
is oblate (Qa(s) < 0) or prolate (Qa(s) > 0). Note in 
particular that the shapes for s = + ^ and for s = +^ 
can be different (see for instance Ref. [ij)- Whenever the 
quadrupole magnetic moment is referred to without ex- 
plicitly mentioning the polarization state, the maximum 
projection is assumed [l|. 



IV. TRANSVERSE DENSITY DEFORMATION 

The interpretation of electromagnetic form factors as 
Fourier transforms of charge and magnetic distribution 
densities is valid only in the nonrelativistic limit. If the 
absorbed photon imparts a significant momentum trans- 
fer to the struck system, the boost of the final state wave 
function relative to the one of the initial state can no 
longer be neglected, which spoils this simple interpre- 
tation. In order to still be able to extract information 
about distribution densities from the measured form fac- 
tors, the concept of a transverse density distribution was 
introduced 041 ■ 

By going to an infinite momentum frame, the depen- 
dence on the longitudinal component of the momentum 
is eliminated, and the deformation is defined in terms 
of densities in the space of the two transverse impact 
parameters, and by. Taking the transverse spin pro- 
jection {s± — ±^,±|) oriented in the x direction, the 
transverse electric quadrupole moment Q^{s±) is (l0| 



Qiis± 



GMi(0)-3eA] 



[Gi;2(0)+3eA] 



}( 



Mi 



(4.1) 



one obtains QX{+^) = -Qil+f )■ 



For Sj_ — ijiit; uuLciiiio i£6A V ' 2^ 

In the previous equation, [Gmi(0) — 3eA] is the electric 
quadrupole moment induced by the moving frame and 
[G£;2(0)+3eA] describes the non-pointlikc structure since 
it is the deviation of Ge2{0) from the "natural value" 
corresponding to a pointlikc particle [lo| . Using the val- 
ues of Eqs. (|l.l[) - (|1.2p . from Eq. (|4.ip we conclude that 
Q^(-|-|) = for a pointlike spin-3/2 particle. 

Similarly, the transverse magnetic octupole moment 



Lattice QCD: 

Quenched \10} 0.83±0.21 

Wilson [10] 0.46±0.35 

Hybrid [10] 0.74±0.68 
Spectator quark models: 

Spectator-S [20| 0.29 
Spectator-SD [22j 0.92 



-3.44 
-3.38 



TABLE I: Transverse electric quadrupole moment Q'^ 



Ok (+ 



2) in units of jfp. 



for the 



becomes [l^l 



+ 1) 



2 {-Gmi(0)-Gb2(0) 
GM3(0) + eA} 



2Ml 



(4.2) 



For s±_ = +i, one has 0^{+^) = -3e'i(+|). A point- 
like spin-3/2 particle has a vanishing transverse magnetic 
octupole moment. 

The authors of Refs. suggest that the electric 
quadrupole and magnetic octupole moments of the A 
should be compared with their natural values, given in 
Eqs. (|l.l|) - (jl.2|) . and that deformation should be defined 
in terms of a positive or negative deviation from those 
reference values. Equations (|4.ip - (|4.2p depend indeed 
on these differences, but what they describe is not the 
full deformation in three-dimensional coordinate space. 
Instead, because the transverse density is defined in the 
two-dimensional [b^^by) plane, they measure an asym- 
metry of the density between the spin direction (along 
the X-axis) and the perpendicular direction (along the 
y-axSs) in the xy-plane. 

To see how much information about the deformation 
of the A is contained in these higher transverse mo- 
ments, we have calculated them for two covariant spec- 
tator quark-diquark models with significantly different 
shapes. The first, model II of Ref. [33|, includes only 
S'-states in the quark-diquark wave function, which is 
therefore spherically symmetric. We call it here model 
"Spectator-S". The second model, presented in Ref. [35j . 
is deformed, because apart from S'-states it includes also 
D-states. We refer to it here as model "Spectator-SD". 

The spherical model Spectator-S yields [201 fo^' the A+ 



G£o(0) = l, Gmi(0) = 3.29, 



(4.3) 



and the quadrupole and octupole moments vanish. For 
model Spectator-SD one obtains [l^] 



Gbo(O) ^ 1, 
G£;2(0) = -1.70, 



Gmi(0) = 3.27 

Ga/3(0) = -1.72. (4.4) 
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From these values we can calculate the transverse elec- 
tric quadrupole and the magnetic octupole moments. 
The results are presented in Table HI together with lat- 
tice QCD data obtained by the MIT-Nicosia group 
with three different methods, for pion masses in the range 
m„ = 350 - 410 MeV. The positive sign of Qi (-h§) 
for all lattice calculations suggests a transverse distribu- 
tion elongated in the spin direction, thus it has a prolate 
(cigar-like) shape. The same deformation is produced by 
model Spectator-SD. 

However, the transverse quadrupole moment for the 
pure iS-wave model Spectator-S is not zero as one might 
expect. Instead, it also predicts a prolate shape, although 
with less deformation. Thus, whereas zero or non-zero 
values of electric quadrupole and magnetic octupole mo- 
ments distinguish clearly between spherical and deformed 
A states, the corresponding transverse moments do not 
provide the same information. This is a quantitative il- 
lustration that the transverse moments are not an unam- 
biguous measure of deformation. 

As for the transverse octupole moment, our result, 
(+f ) < 0: suggests a deformation perpendicular to 
the spin axis, thus an oblate (pancake-like) shape. This 
is true for both Spectator quark models, with and with- 
out ZJ-states. The numerical values arc very close, which 
means that docs not discriminate much between 
models with spherical or deformed wave functions. 

At the moment, no lattice calculations of are avail- 
able. But we can use the form factor data of the Ade- 
laide group at Q2 = 0.23 GcV^ in Eq. g^), replac- 
ing CA by Geo{Q^), for a rough estimate. We obtain 
CaC+I ) = (-23.8 ± 22.3) 2]gT, which is consistent with 
our result — although with large statistical uncertainty — 
and also suggests an oblate shape. For a more rigorous 
comparison with our predictions we have to wait for fu- 
ture lattice calculations. 



DEFORMATION IN SPECTATOR QUARK 
MODELS 



In this section, we calculate the A charge densities 
from wave functions obtained in a covariant spectator 
quark model, and we illustrate to what extent the distor- 
tion caused by the D-wave contributions manifests itself 
both in the momentum-space and coordinate-space den- 
sities. There is no need to make use of electromagnetic 
form factors to characterize deformation in this case, be- 
cause the densities are calculated directly from the wave 
functions. Thereby we sidestep the usual problems in 
relating densities and form factors, namely that the lat- 
ter involve wave functions in different reference frames. 
The former is the Fourier transform of the latter only in 
the nonrelativistic limit, whereas we are interested in the 
general, relativistic case. 

Because the covariant spectator theory is more nat- 
urally formulated in momentum space, in the follow- 
ing we start our discussion of A charge densities in the 



momentum-space representation where the wave func- 
tions were developed [l2|, [s^l • Then we perform a Fourier 
transform of the wave function and discuss the densities 
in coordinate space. 



A. Momentum space 

The components of the A wave function of Eq. (|2.ip 
can be written [s^l as 

«'s(P, k- s) = -Vs(P, fc)(£p)" Wa(P, s) , (5.1) 
*I33(^, fc; s) = ^D3(P, k)^D3{P, k; s) , (5.2) 
^Di{P,k;s)=i'Di{P,k)'S>Di{P,k;s), (5.3) 

where the isospin state, which is a common factor in all 
wave functions, is omitted, and we use the notation 

'^>D(2S)iP,k;s) ^ -3{e*pr irs)^f;V^''u,iP,s), (5.4) 

where S = i , | is the core spin and P the four- 
momentum of the A. In Eqs. (|5.ip and (|5.4p . Sp is the 
diquark polarization vector in the fixed-axis representa- 
tion [2^ . In the last equation, Vs is a projector onto the 
|, and V is the spectator Z?-state 



state 5' = i or S" 

operator. More details can be found in Ref. [34 

The scalar functions ipSi i^oa, and ipoi regulate 
the momentum distribution of the quark-diquark sys- 
tem. Because both the A baryon and its diquark con- 
stituent are on-mass-shell, these functions depend only 
on (P — fc)^. It is convenient to express them in terms of 
the dimensionless variable 



X ■ 



(Ma - - (P - fc)^ 
MattIs 



(5.5) 



where rus is the diquark mass. Following Ref. |35| . we 
use the paramctrizations 



^s{P,k) = 

i^D3iP,k) 

i^Di{P,k) 



Ns 



_ Nd3 

ndi r 1 



(5.6) 
(5.7) 
(5.8) 



The momentum range parameters ai {i = 1,...,4), 
given in units of mg, determine the long- and short- 
range dependence of the wave functions in coordinate 
space. The normalization constants Ns, N03, and Noi 
arc determined by the conditions J^, |?/;5(P, fc)p = 1, 
/^k4|V'B3(P,fc)|2 = 1, and J^k^\i,Di{P,k)\^ = 1, where 
P = (Ma, 0) is the total four- momentum of the A baryon 
in its rest frame. The coefficient X^i is determined 
through the orthogonality between the A and nucleon 

states [H,!!!!!!. 

In this work, we use the model of Ref. [1^, with the 
parameters ai = 0.3366, 02 = 0.35054, 0:3 = 0.33773, 





1 2 3 4 5 

k (GeV) 



FIG. 2: Comparison of the three contributions to the total 
momentum-space density PA(k, s) in Eq. 115.1311 in units of GeV~^. 
The solid line represents the symmetric contribution, p£(k), the 
dashed and dotted lines show the coefficients of Y2o{z) propor- 
tional to ^s{P}k)ipD3{P}^) and \iiDl{Py ^)\'^ , respectively. In all 
cases, the common factor ba = 1, and only the absolute values arc 
plotted. 



and a4 = 0.34217. One obtains Am = 1.031898, and 
the Z3-state admixture coefficients are a = 0.08556 and 
b = 0.08572. 

Similarly to the case of the nucleon in Ref. [l^, the 
momentum-space charge density of the A in its rest frame 
is defined as 

PA(k,s) =^*^(P,fc;s)j,*A(P,fc;s), (5.9) 

As 

where s = ±i, ±| is the spin projection of the A state, 
and jq = 3ji = ^ + the charge operator [l^, [33| . 
Remember that implicitly the wave function '^a{P, k; s) 
depends also on the diquark polarization A,,. 



Substituting (j2.ip into (|5.9|) one gets 
PA(k,s) = ArVA,s(k,s) 

+2aN^PA,SD3{Ks), (5.10) 

where 

PA,xiKs)=J2^UP^k;s)jq^x{P,k;s), (5.11) 

As 

for X = S,D3,D1, and 

PA.SD3 (k, s) - ^ *t^3(F, k; s)jq^s{P, k; s). (5.12) 

As 

The S and D3 states can either be in the initial or final 
state, hence the factor of 2 in the last term of Eq. (|5.10p . 

After performing the spin and isospin algebra, it is 
helpful to isolate the spherically symmetric contribution 
from the angle-dependent terms in the total density, 
and one obtains 

PA(k,s) = pi(k) 

+2eAaNy,is)k' [^siP,k)^D3iP,k)] ^20(2) 
-eAb^N'f,{s)k^\i^Di{P, k)\%o{z), (5.13) 

where z = cos 9, the function Y2o{z) = ^(3z^ — 1) is 
proportional to the spherical harmonic Y20, and 

pi(k) = caN^ X (5.14) 
[\MP,k)\^ + a'k''\^D3iP,k)\' + b^k^\^Di{P,k)\^] . 

In the rest frame, tpx depends only on k^. The second 
term in (|5.13p is the overlap between the 5" and the D3 
state. It does not vanish because the states have the 
same core spin-3/2. The Dl state on the other hand 
is always orthogonal to S and D3 because its definition 
includes a spin-1/2 projector. The last term comes from 
the overlap between the initial and final Dl states. The 
two last terms of Eq. ()5.13p vanish when the angular 
integration is performed. 
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FIG. 3: Contour plots of momentum-space charge densities of 
the A in the k^-kz plane, in units of GeV~^. The top panel shows 
the total density pA (k, + §)• The bottom panel isolates the angle- 
dependent part p^{k) induced by the D-states. Note the difference 
of the density scales in the two panels. 

The radially symmetric part of the density can also be 
written as 

pi(k) = i [pa (k, +1) + PA (k, +i)] , (5.15) 
and the angle-dependent asymmetric component as 
pKk) = i[pA(k,+|)-pA(k,+i)] 

^b^eAN^k^\ijDi{P,k)\^Y2o(z). (5.16) 

Using (|5.14p and (|5.16p we can rewrite pA(k, s) as 

PA{k,s)=pi{k) + Ms)pi{k). (5.17) 

From this equation and from fs (±^) = —fs (if) it is 
clear that the deformation density for s = +i has always 
the opposite sign of the one for s = +|. 




-OA -02 0.0 02 04 

ft, (GeV) 




-6.4 -6.2 6.0 o!2 04 
iGeV) 

FIG. 4: Contour plots of momentum space charge densities of the 
A in the kx-k~ plane, in units of GcV~'^. The upper panel shows 
the Dl — Dl part, the lower panel the S — D3 part of p^ (k). Note 
the difference of the density scales in the two panels. 

For small D-state admixture coefficients a and 6, the 
S to D3 transition term dominates the asymmetry. This 
is consistent with a calculation of the form factors in first 
order in a and b where the D3 state is responsible for a 
non-zero electric quadrupole form factor (2ll. [22j|. 

The factors tpsi'Ds and tpf^-^ in Eq. (|5.13p are always 
greater than or equal to zero, and therefore the D3- and 
Z?l-state contributions to the deformation enter with op- 
posite signs. The overall factor multiplying Y2o{z) can 
have either sign, depending on the specific parametriza- 
tion of the wave functions and on the spin projection s. 

To illustrate the deformation of the A graphically, in 
Fig. [T] we show a polar representation of pA(k, s), for 
k in the kr^-k^ plane (the densities are invariant under 
rotations about the fc^-axis). The positive k^ direction, 
corresponding to polar angle 6 — 0, points upwards. For 
fixed values of k = |k|, the length of a straight line from 
the origin to a given point on a displayed curve is the re- 
spective density, and its angle with the upward direction 
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is the polar angle 0. 

In this representation, a distribution p^(k,6]s) with 
no ^-dependence yields a perfect circle. If pA{^,d',s) is 
larger in the than in the direction, the deformation 
is oblate, whereas a larger p/^(k,9;s) in the k^ than in 
the kx direction corresponds to a prolate deformation in 
momentum space. 

The three panels of Fig. [l]show pA(k, s) for k = 0.2, 
1.0, and 2.0 GeV, respectively. In each case, the dashed 
line represents s = +|, and the dotted line s = 
The perfect circle (solid line) shows the symmetric dis- 
tribution p^Ck). In the case of the s = -|-| density, the 
deformation is prolate for the smaller momenta k=0.2 
and k=1.0 GeV, but is oblate for k=2.0 GeV. Recall that 
the deviations from a spherically symmetric density for 
s = +^ and s = -|-| are equal but with opposite signs. 

This change in the shape of the deformation with in- 
creasing momentum can be understood from the behav- 
ior of the S — D3 and Dl — Dl terms in the density 
of Eq. (|5.13|) . Figure [2] shows the k-dependence of the 
magnitudes of 2aN'^]s.^ipsipD3 and b'^k'^N'^ipj^^, the fac- 
tors that multiply Y20 with opposite signs, in compari- 
son with the symmetric term /9^(k). For small k, the 
S — D3 term dominates strongly over the Dl — Dl term, 
which is expected already because it is of first order in the 
small D-wave admixture parameter a, whereas Dl — Dl 
is quadratic in b. As k increases, the S — Di contribu- 
tion falls faster than Dl — Dl, due to the faster falloff 
of the ^-state wave function, and becomes equal in mag- 
nitude at around k ~ 1.2 GeV, where they cancel. For 
larger values of k Dl — Dl dominates over S — D3 and 
changes the shape from prolate to oblate. However, the 
coefficients of Y20 at such high momenta are already very 
small (note the logarithmic scale in Fig. [5]). 

Figured] shows also that, at high momenta, the spheri- 
cally symmetric part of the density, p^ (k) , is itself domi- 
nated by its _D1 — _D1 component. Comparing Eqs. ()5.13|) 
and (|5.14|) we sec that the magnitude of the Dl — Dl co- 
efficient of Y2o{z) in (|5.13|) . represented by the dotted 
line in Fig. [2l is the same as the Dl — Dl contribution to 
the spherically symmetric part of the density in (|5.14p . 
With increasing momentum k, the dotted line seems to 
converge to the solid line of the total symmetric contri- 
bution. However, it reaches only about 90% of p£(k). 
the small D3 — D3 contribution being responsible for the 
remainder. This ratio of the Dl — Dl contribution to the 
total spherically symmetric part can be obtained using 



the asymptotic ratio — . — ^. 

Vol 1 — ADl Ndi 3 

Figure |3] shows contour plots of momentum space 
charge densities in the k^-kz plane. The deformation 
of the total density pA (k, -|-|), displayed in the upper 
panel, is barely visible in this plot, because the spher- 
ically symmetric contribution dominates strongly. The 
lower panel shows only the much smaller asymmetric part 
p^(k), which enhances the density along the /cz-direction. 

The origin of this enhancement is analyzed in Fig. |4l 
where the upper panel shows the Dl — Dl, and the lower 



panel the S — D3 contribution. The Dl — Dl channel 
density is oblate, but at the low momenta shown in the 
figure it is overwhelmed by the more than one order of 
magnitude larger prolate S — D3 density deformation. 

In summary, at low and intermediate momenta the mo- 
mentum space density distribution shows a small prolate 
deformation for the A state with spin projection s ~ 
and a small oblate deformation for the s = +^ state. For 
high momenta the shape of the deformation is reversed, 
but this is hardly relevant because the density is already 
negligibly small in this region. 



B. Coordinate space 

We proceed now to the calculation of charge densi- 
ties of the A in coordinate space. All densities will be 
determined in the rest frame of the A, where the total 
four- momentum is P = (Ma,0). Using the modified 
momentum-space wave function 



(5.18) 



where the relativistic phase-space factor is absorbed into 
the definition of the wave function, we can write the 
Fourier transform in the same form as for nonrclativistic 
wave functions: 



^A(r;s) = 



(2^ 



e'^''^Aik;s). 



(5.19) 



With this convention for the factors of 27r, the inverse 
transform is 



*A(k;s) - Jd^r e-''' '-^'A(r;s). 



(5.20) 



In complete analogy to the momentum-space expression 
(|5.9p . the charge density in coordinate space is given by 



PA{r; s)=Y^ *i(r, .s)jg*A(r, s). 



(5.21) 



As before, one can decompose ^'A(r;s) into angu- 
lar momentum components S, DS, Dl. In the fixed- 
axis polarization state basis we are using, there is no 
fc-dependence in the diquark polarization vector (e^)^. 
This makes the Fourier transform of the S'-state (|5.ip in 
the rest frame particularly simple. 



^s{r; s) = -i?s(r)£>"(^, s). 



(5.22) 



where we have introduced the shorthand e* to represent 
(e%) _ , and 



d^k 
(2^ 



-V's(P,fc) 



(5.23) 



with r ~ |r|. The factor e* depends on the diquark po- 
larization As, which is not shown explicitly. The isospin 



states are not affected by the transformation and are also 
suppressed for simplicity. There is no angle-dependence 
in the rest frame wave function ■(/'^(P, fc), therefore Rs 
is also an S-wave, depending only on r. Boosting the 
wave function to another frame would induce an angle- 
dependence into tpsiPj k) and consequently also into Rs- 
However, this kind of relativistic deformation due to 
Lorentz contraction of the system along the direction of 
motion is a separate issue. We are interested in intrinsic 
deformations of the A, which are already present in the 
rest frame. 

Using the familiar expansion of a plane wave into par- 
tial waves, 

-t-oo I 

e^k-r^4^^ ^ fyUmrMjii^^), (5.24) 

1=0 m=-l 



where ji are the spherical Bessel functions, Eq. (|5.23p 
becomes 



(5.25) 



The general structure of the complete Z?-state wave 
functions (without isospin) can be written as 



*i)(k;s) 



<pDik,s) 



(5.26) 



where D stands for Dl or D3. $D(fc, s) is the spin wave 
function of the _D-statcs. In the rest frame, it can be 
represented as [131 

$Z3l(fc,s) = +V^khl 

X ^(2mz;isi||s)r2„,(k)t/"(P,si) 

m,si 

X J2 (2m,;|si|fs)y2™<(k)^^"(P,si) 



(5.27) 



where U°' is the spin- 1/2 state |33l. 134 



C/"(P, s) ^ -^75 (7" - ^) uiP, s). (5.28) 

The angle dependence is contained exclusively in the fac- 
tors 12771, (k), and the spin states are completely indepen- 
dent of k. 

The 7D-state wave functions in coordinate space are 
then given by the Fourier transform, which yields 



^'m(r;s) = +V47r 

X ^(2m,;isi|f,s)3^-|(r) [slU"{P,s,)] 



*D3(r;s) 



mjsi 



mjsi 

(5.29) 
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FIG. 5: Radial S- and Z?-state wave functions of the A in coor- 
dinate space, calculated through Fourier transforms according to 
Eqs. Is:^ and ll532t . 



where 



(27r)3 



(5.30) 



Again, in the rest frame ipD{P,k) is also angle indepen- 
dent, and the Fourier transform simplifies to 



-i?Z5(r)r;™,(f), 



(5.31) 



where 



RDiT)=4n 



(2^ 



, 2 . /, ■, i'D{P,k) , . 
k j2(kr) — 1=^, (5.32) 



and the minus sign factored out in (j5.3ip comes from the 
in dEHl). 

We calculated the functions Rs, Rds, and Rdi numer- 
ically, and the results are presented in Fig. [5] The S-state 
dominates at small distances, but the two D-state wave 
functions become comparable in size around r w 1 fm 
and dominate for larger values of r. The D-waves start 
out with opposite signs, but the Dl wave changes sign at 
r « 1.4 fm. 

As r goes to zero, the relativistic I?-state wave func- 
tions are weakly singular, namely i?z)(r) cx r~^/^. This 
behavior, which is reminiscent of the singular radial de- 
pendence of the Dirac wave functions of the hydrogen 
atom 41[ , does not cause any problems because the den- 
sities remain integrable. The origin of these singularities 
is the slower falloff with increasing relative momentum 
of the relativistic momentum-space wave functions com- 
pared to nonrelativistic wave functions. If we calculate 
the Fourier transform of the D-state wave functions in 
the nonrelativistic limit, at small r we obtain the regular 
behavior i?_D(r) cx r". 

The total coordinate-space charge density is 



PA(r,s) ^ iV^PA,s(r,.s) 

i^N^PA,D3{^, s) + b^N^PA,Diir, s) 



+2aN^PA,sD3{Ks). 



(5.33) 
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FIG. 7: Comparison of the three contributions to the total 
coordinate-space density p/\(r,s) in Eq. 115.461 1 in units of GeV^. 
The sohd line represents the symmetric contribution, p£(r), the 
dashed and dotted lines show the coefficients of Y2o{z) propor- 
tional to RsRd3 and R^i, respectively. In all cases, the common 
factor e/^ = 1, and only the absolute values are plotted. 



The various components are defined in analogy with 
([5lT|) and dSH]). 

The density associated with the iS-state is 



PA,s{r,s) = ear's ■ 

For the other cases, one uses 



(5.34) 



(5.35) 



M„(P, s)u"(P, s) = L7„(P, s)C/"(P, s) = -1 , (5.36) 
UaiP,s)U°'{P,s) = U°'{P,s)ua{P,s) =0, (5.37) 



assuming the same polarization in the initial and final 
states. The results for D3, Dl, and the transition 5 to 



for s 



-|, and 



D3 I 
2 

Dl 



-eARDs.RsY2o{z), 



PA,_Di(r, + 
PA.SDsfr, 



caR 



Dl 



= eARD3RsY2aiz) ■ 



(5.38) 
(5.39) 
(5.40) 

(5.41) 
(5.42) 
(5.43) 



for s — The function Y2o{z), with z = cos 6, was de- 
fined previously, but the angle 9 is now to be understood 
as the angle between f and the z-axis. More details are 
given in Appendix Rl 

The total density becomes 



PA{r,s)=pi{r) + fs{s)p^{r), 



(5.44) 



where 



pi(r) = N^CA [Rl + a^Rls + b^Rli] (5.45) 

is again the angle- and spin-projcction-indcpcndent con- 
tribution, and 



'2eAaN^RsRD3Y2oiz) 
-eAb^N^Rl^Y2o{z), 



(5.46) 



the angle-dependent asymmetric component. Equa- 
tion (|5.46p is the coordinate-space analogue of Eq. (|5.16p . 
The different sign of the term containing the S and D3 
wave functions is due to Eq. (|5.3ip where i?_D(r) is defined 
such that it does not contain the factor iK 

Again there are two independent terms proportional 
to Y2o{z) that cause deformation, one associated with 
the Dl state, and another with a S* to D3 transition. 
The main difference to the analogous momentum-space 
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FIG. 8: Contour plots of coordinate-space charge densities of the 
A in the x-z plane in units of GeV^ . The top panel shows the total 
density (r, -|-|) . The bottom panel isolates the angle-dependent 
part (r) induced by the D-states. 



expression is that both terms come now with a minus 
sign. Because Rdi enters squared, the total effect of the 
two iJ-state contributions can be inferred from the sign 
of the wave functions Rs and Rds and the sign of a. 

When a is positive, as in model Spectator-SD, and Rs, 
Rd3 have the same sign, the two terms in Eq. (|5.46|) have 
also the same sign and reinforce each other's contribution 
to the deformation. This is the case in the region up to 
r = 4 fm shown in Fig. [5l where both Rs and Rds sltc 
positive. We checked that both Rs and Rd3 stay positive 
in that region when the wave function parameters are 
varied within a broad range, such that the direction of 
the deformation remains a robust result. 

For s = +1 this means that Y2o{z) is multiplied by 
an overall negative factor, implying an oblate shape, as 
shown in Fig.[6]for r = 0.2, 1.0, and 2.0 fm. For s = -f i 

the coefficient of Y2a{z) has the opposite sign, and the 
deformation is prolate. At r = 0.2 fm the deformation is 
too small to be visible in the graph, but it becomes more 



pronounced as r increases. 

Figure [7] shows the magnitudes of the individual co- 
efficients of Y2o{z), together with the radially symmetric 
y5^(r). The S* — D3 contribution dominates over Dl — Dl 
at small r, because it is of first order in the small Z?-state 
admixture coefficient a, whereas Dl ~ Dl is oi second 
order in b. However, with increasing r the S'-wave falls 
more rapidly than the D-waves, and the £)l-term domi- 
nates for r > 3 fm. 

Finally, Fig.[8]presents contour plots of the charge den- 
sities of the s — -|-| state in the x-z plane. The upper 
panel shows the total density, the lower panel only the 
much smaller asymmetric contribution which is respon- 
sible for the oblate shape. 

Our results for the shape of the are in agreement 
with those of pre vious studies, such as the results re- 
ported in Ref. [IJ] using the cloudy bag model^nd the 
constituent quark model calculations of Ref. [l|. The 
oblate shape for A+ is also consistent with the negative 
sign of Ge2{0) found in lattice QCD simulations [loll40|. 

The formalism presented here can be applied to other 
systems. For instance, from an analysis of the form fac- 
tors of the r2~ baryon, it can be concluded that the charge 
density distribution of the £7" also has an oblate shape 

M. 



VI. SUMMARY AND CONCLUSIONS 

The A is the lowest-mass baryon that can possess a 
non- vanishing electric quadrupole moment. In a nonrel- 
ativistic framework, the electric quadrupole moment can 
be used as an indicator for a particle's deviation from 
a spherically symmetric shape. However, in the general, 
relativistic case the connection between shape and higher 
spectroscopic moments is more complicated, which led to 
the proposal of alternative methods to measure deforma- 
tion. 

One of these methods suggests to extract information 
about deformation from transverse densities, calculated 
in the transverse impact-parameter space (hx,by) in the 
infinite-momentum frame. It has the advantage that the 
transverse density moments and are zero for 
pointlike particles, and therefore can be used to measure 
the extension of the particles in the impact-parameter 
space. However, they do not allow the classification of 
the shape in the particle's rest frame. 

In this work wc used a different relativistic formalism, 
the covariant spectator theory, to investigate the rela- 
tion between moments of charge or magnetic density dis- 
tributions and the intrinsic shape of these distributions 
in the A baryon's rest frame. We used two covariant 
quark-diquark momentum-space wave functions for the 
A, one consisting of pure 5- waves only, called "Spectator- 
s'' , and another which includes S- and _D-waves, called 
"Spectator-SD" . The electric and magnetic moments and 
form factors can be calculated directly from these wave 
functions, as well as the momentum space densities. The 
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coordinate space densities are then obtained from the 
Fourier-transformed wave functions. 

We arrived at the following results: 

For the S'-wave model Spectator-S one obtains 
G_E2(0) = Ga/3(0) ~ 0, and the electric and magnetic 
density distributions are spherically symmetric. 

The D-wave admixture in model Spectator-SD on the 
other hand produces a spatial deformation of the A den- 
sity distribution, and the quadrupole and octupole mo- 
ments become Ge2{0) = —1-70 and G'a/3(0) = —1.72, re- 
spectively. The negative value of the quadrupole moment 
corresponds to an oblate density distribution in coordi- 
nate space (prolate in momentum space). We conclude 
therefore that the higher moments are good indicators 
that allow to distinguish deformed from spherically sym- 
metric systems. 

Using the same wave functions and their respective 
electric and magnetic moments, we also calculated the 
corresponding transverse density quadrupole and oc- 
tupole moments (Q^ and O^). For each moment, the 
obtained values (sec Table |I| for models Spectator-S and 
Spectator-SD are non-vanishing and of the same sign, 
and thus do not show a clear distinction between spheri- 
cally symmetric and deformed cases. Since the moments 
are non-zero one can conclude that the system is not 
pointlike in impact parameter space, but it is not clear 
how further information on its shape could be extracted. 
In this sense, the usual three-dimensional density dis- 
tribution complements the information contained in the 
transverse densities, and it is also closer to our intuitive 
notion of deformation. 

For the specific case of the baryon with spin projec- 
tion s = -|-|, the covariant model Spectator-SD predicts 
an oblate shape of its density distribution in coordinate 
space. This is in agreement with previously obtained 
results, both from other quark model calculations and 
lattice QCD simulations. 

All model parameters were determined through fits to 
the available lattice QCD data for the 77V — > A transition 
form factors at large pion mass, where the uncertain pion 
cloud effects are minimal (the experimental data are also 
well predicted) (ssj . In particular, the coefficients a and 
b are determined by the lattice transition form factors 
for ~ 0, namely a is determined by by the Coulomb 
transition form factor G^{Q'^) and b by the Coulomb 
transition form factor Gq^Q"^). Further improvements in 
the statistical quality of the lattice data might alter the 
magnitudes of a and b and therefore the extent of the 
deformation we predict for the A, but the signs are not 
in doubt. In this sense, the nature of the deformation, 
namely that it is oblate rather than prolate, is a robust 
prediction of our model. 
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Appendix A: Explicit expressions for 4'i3i(r;s) and 
*i33(r;s) 

We list here more explicit expressions of the 
coordinate-space A wave functions 'I'£)i(r;s) and 
\I'£)3(r; s). Performing the sums in Eq. (|5.29p we get 



^Di (r;+l] 



Wtt. 




and 



(r;+|) 



^yil{r)el [/" (P,-i),(Al) 



3^il(r)4 [/" (P,-i).(A2) 

-3^&3(r)£>" {P,+l) 
-3^+i(r)e>" (P,+i) 
'.y+l{v)elu-{p,-l), (A3) 



/ AlT - 

-\lfyilir)e:u'^{P,-l)- (A4) 

In the calculation of the density distributions from the 
coordinate wave functions we used the following rela- 
tions: 

4^ 1^ \Y,+,ir)f + ^ \Y2+^{rf + J \Y,oir)f^ = 1 

4^ 1^ \Y2+2{r)\' + I \Y,oir)\' + ^ \Y,^,{r)f^ = 1, 

(A5) 

4^ 1^ \Y,+,ir)\' + 1 |r2+i(f)l'} = 1 - ^ (3^' - 1) 

4^ 1^ |l2+i(f)|^ + ^ |r2o(f)|'} = 1 + 1 {3z' - 1) . 

(A6) 
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